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Abstract 

It is shown that the inclusive rate of the rare weak radiative decays B — * X s 7 
contains a series of nonperturbative corrections, whose 'short distance' scale is set by 
m" 1 , rather than by m^ 1 . The first correction in this series is expressed through 
the chromomagnetic interaction of the b quark inside the B meson and the relative 
magnitude of the effect is determined by the ratio (B\ba • Gb\B) jrr? c . Though the 
magnitude of this first correction is suppressed by a numerical coefficient, the sensitivity 
of the decay rate to the distance scale m~ l may significantly limit the accuracy of purely 
perturbative predictions for the rate. 



The rare radiative decays B — > X s 7 associated with the underlying FCNC quark process 
b — > S7 are well known to be sensitive to the top quark mass as well as to a possible 'new 
physics' that might exist beyond the Standard Model. The current state of the experimental 
knowledge about these processes and of their understanding within the Standard Model can 
be summarized by the measured valuei of the inclusive branching ratio of these decays 

B(B -> X s 7) = (2.32 ± 0.57 ± 0.35) x 10" 4 (1) 

and by the perturbative Standard Model quark decay rate, whose calculation has been 
recently taken @ to the next-to-leading order (NLO) in a s ln(mw/j7ij,), 

B(b^X sl ) = (3.28 ±0.33) x 10~ 4 . (2) 

Additionally, it has been arguedB that the nonperturbative effects, which distinguish be- 
tween the heavy quark decay and the decay of a heavy meson, are of relative magnitude 
A.Q CD /ml. Thus in the subsequent studies (see Ref. and references therein) it has been 
assumed that the perturbative total quark decay rate approximates the total inclusive meson 
decay rate with accuracy better than 10%, though the 0(l/m%) nonperturbative effects are 
naturally expected to be significant in the photon spectrum near the endpointS. 

The purpose of this paper is to point out that there in fact exists a series of nonper- 
turbative corrections to the inclusive rate of the decays B — > X s 7, whose 'short distance' 
scale is set by m" 1 , rather than by . This series thus generates an expansion for the 
corrections to the decay rate in powers of AQCD/ m c and of AQCD m b/ m c- The first term in 
this expansion is 0(Aq CD /m 2 l ) and is found in terms of the known chromomagnetic energy of 
the b quark. The subsequent terms, however are expressed trough unknown matrix elements 
of operators of higher dimension. The relative magnitude of the calculable first correction 
to the total decay rate is found here in the lowest order in a s and is given by 



6T(B^X sl ) 1 » 



-0.025 , (3) 



T(b^X sl ) 27 C 7 m 2 c 

where \j? is the standard parameter || for the strength of the chromomagnetic interaction 
of the b quark inside the hadron: ^ = (B\bga^G% (A a /2) b\B)/2 = 3 (M|, - M|)/4 « 
0.4 GeV 2 . (The nonrelativistic normalization for the heavy quark states is assumed through- 
out this paper, so that (B\tfb\B) = 1.) The quantity C7 in the denominator in eq.([|) is the 
coefficient of the operator 

g 

p 7 = 7F~9 rn b (s L b R ) (4) 
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in the perturbative effective Lagrangian for the b quark decay normalized at a scale /i ~ m b : 

L eff = ^V: s V tb jzC^)P^) , (5) 
v 2 i=l 

and the values C-j = —0.30@ and m c « lAGeV are used in the numerical estimate in 
eq.@Q. Clearly, the results of the present paper, in particular eq.(|3]), are at variance with the 
generally accepted viewpoint that the relative magnitude of the nonperturbative corrections 
to the inclusive decay rate of B — ► X s ^ goes to zero in the limit mj, — > oo. The latter 
behavior would be true if the relevant contributions to the B meson decay were exhausted 
by the effective Lagrangian in eq.(|5J), which is the case for the perturbative b quark decay. 
At the level of nonperturbative effects there are additional terms in the effective Lagrangian, 
among which there are the ones of the type considered in this paper, whose distance scale 
is set by m' 1 , so that their (relative) contribution does not vanish in the limit m b — > oo. 
The correction in eq.(H) is small due to a small numerical factor. Apriori there is no reason 
to expect that similar coefficients are small in subsequent terms. Moreover, the subsequent 
terms contain average values of operators of the generic form (b(q D) n Gb) / 'm 2 c n+2 , where q 
is the momentum of the gluon and the covariant derivative D is acting on the gluon field 
tensor G. After averaging over the phase space in the decay these contributions produce 
series in (Ag CD /ra;!) (AQ C D m b/ m l) n - Since in reality the parameter AQ C D m b/ m l is of order 
one, these terms bring a considerable uncertainty in predictions for the decay rate. 

The consideration of the nonperturbative effects in the inclusive decay rate is performed 
using the standard method© 0. Namely, both mj and m c are assumed to be large: m&, m c ^> 
Aqcd, and one performs the operator product expansion in the inverse powers of the heavy 
quark masses for the effective operator 



21m 



d'xe^T{Lt ff (x),L eff (0)} 



(6) 



The total decay rate for a heavy hadron Xj, is then given by the diagonal matrix element of 
T over X b : 

T(X b ) = (X b \T\X b ) . (7) 



1 The notation and conventions of Ref.[g[ are used here as well as the numerical value of C-r(mb). Fur- 
thermore, eq.([|) is written in the approximation where T(b — > X s j) — T(b — > sj) and the latter rate is 
generated only by the operator P-j. In NLO the rate also receives a small contribution from other operators 
in L e ff. However, since the correction is calculated here only in the leading order in a s , taking into account 
those additional contributions to T(b — * X s 7) would exceed the accuracy of the present calculation. 
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Different types of decay are separated within this method by picking up the appropriate terms 
in L e ff to correlate in eq.(|^). The leading term in the expansion generates the perturbative 
'parton' decay rate, while the terms suppressed by inverse powers of heavy quark masses 
describe the nonperturbative corrections to the decay rate. 

In the leading logarithm approximation (LLA) the perturbative quark decay rate of 
B — > X s 7 is obtained by retaining in L e ff only the term with the operator P-j in eq.@ and 
substituting it in the correlator T in eq.(fD with the result 

T s? = ^ G F m l\ V t s V tb \ 2 \C 7 \ 2 (bb) . (8) 

This expression gives the standard perturbation theory formula for the inclusive rate of 
B — > X s 7 upon using noticing that up to terms 0(m^ 2 ) one has (B\(bb)\B) = 10. 

The reasoning here for terms suppressed by only m~ 2 rather than by m^ 2 is based on 
considering the 'gluon - photon penguin' type of mechanism, which provides additional to 
eq.(|5]) terms in the effective Lagrangian. This mechanism arises through graphs of the type 
shown in Fig. 1. The gluon momentum k in this graph is assumed to be small, and the 
expansion in this momentum gives rise to a series of operators in L e ff with the gluon field 
strength tensor G, and its derivatives. One can readily see, that already starting from first 
such term the loop integration for the coefficients in this expansion is convergent in the local 
four-fermion limit of the weak interaction. In particular, the corresponding graph with the 
top quark becomes irrelevant. In the local limit of the weak interaction the loop with the 
c quark is kinematically equivalent (after the Fierz transform of the four-fermion vertex) 
to the c quark loop for coupling of axial current to two vector currents. Since the photon 
momentum q is on-shell: q 2 = 0, and the gluon momentum is infinitesimal, the external 
kinematical invariants for the loop, q 2 , k 2 , and (q + k) 2 are either zero or infinitesimal. Thus 
the coefficients of expansion of the loop in powers of k are determined by the mass of the 
quark in the loop, i.e. by m c . Furthermore, an expansion of the loop in the invariant (q + k) 2 
gives rise to terms with powers of the ratio (q ■ k)/m 2 , resulting in the operator terms of the 
form (b (q D) n G b)/m 2n+2 , which give rise to corrections to the rate of the relative magnitude 
^qcdIO (A Q CDm b /m 2 c ) n . 

By calculating the graph of Fig. 1 and the one with the photon and the gluon permuted, 
and expanding to the first order in the gluon momentum k, one explicitly finds the first term 
in this series of corrections as an additional contribution to L e ff of the form 
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where Q c = 2/3 is the electric charge of the charmed quark and the expression is written 
under the convention about the sign of Gp as in Ref. 0, i.e. where the bare four-fermion 
Lagrangian is written as L(b — > ccs) = -(AG F /V2) V* s V cb (s L 7 M c L ) (c L 7 M 6 L )- 

When correlated according to eq.(|6|) through the on-shell (massless) s quark and the 
photon with the term with the operator P-j in L e ff of eq.([|) the newly found term L^j^ 
gives an interference contribution to the effective operator T. This contribution is found 
here as 

Tff = Gl m\ lite [V* V cb V ts V* C 7 ] ^fjj® ■ (10) 

Ignoring the terms in the unitarity relation for quark mixing associated with the uu pair 
contribution, which are of the relative order (| V^/V^l) sin 9 C , one can write V* s V cb = —V t * V t b- 
Then a comparison of the equations (|T0|) and @ gives the expression in eq.(§) for the relative 
nonperturbative correction. 

Proceeding to a discussion, we first remark that the ignored contribution of the uu pair 
in fact leads to a hopefully small, but presently uncalculable effect in the rate of B — > X s 7. 
The smallness is obviously related to the suppression in the mixing, while the uncertainty of 
the contribution arises, when one considers the same loop as in Fig. 1, but with the u quark 
instead of the charmed one. In no approximation can the mass of the u quark be considered 
as large. Thus this mechanism gives rise to an essentially long distance contribution to the 
inclusive rateQ. 

Returning to the discussion of the corrections arising from distances C^m" 1 ), it can be 
noticed that in general the effect of such distances can also significantly affect the spectrum of 
inclusive photons through enlarged distance scale for the amplitude of the photon emission. 
However any detailed consideration of this effect is beyond the scope of the present paper. 

The particular numerical value in eq.(§) for the relative correction to the inclusive rate 
is most certain to be modified by the subsequent terms with derivatives of the gluonic 
tensor. Apriori those terms have no parametric smallness in comparison with the expression 
in eq.(|3|), since the additional factors in those terms are given by powers of the parameter 
Aqcd mb/ml ~ 0(1). The principal difficulty of evaluating those subsequent terms is in that 
we have no knowledge of the matrix elements of operators with derivatives of the gluonic 
tensor over the B meson. Therefore this series of the corrections puts an essential limit on 

2 Possible long distance effects were in fact discussed within a different technique in the literature^ . 1 
believe that the used here systematic method of OPE provides a more reliable approach to calculable terms 
and leaves undetermined only a small contribution due to the light uu pair, which is strongly suppressed by 
the relative factor (IK^/V^I) sin6> c « 0.02. 
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the present calculability of the decay rate. 

It can be noted that the presence of a contribution in the rate due to the chromomagnetic 
interaction of the b quark can in principle be verified experimentally by measuring also 
analogous inclusive rate for the A& baryon: r(A& — ► X s j). In A& the spin of the light 
component of the baryon is zero, thus the average of the chromomagnetic operator vanishes. 
Therefore for A& the discussed correction is absent, and the difference between the inclusive 
decay rates T(B — > X s 7) — r(A& — > X s 7) is dominated by the spin- dependent terms in the 
discussed corrections with further small corrections being 0(m^ 2 ) and analogous to those 
found@ for the dominant decays. Thus this difference in rates can serve as a measure of 
importance of the nonperturbative corrections coming from distances (^(m^ 1 ). 

In summary. It has been shown that there exists a set of nonperturbative corrections 
to the inclusive rate of the decays B — ► X s 7, which are determined by the distance scale 
(^(m" 1 ). The relative magnitude of these corrections is not vanishing in the limit m;, — > 00, 
provided that m c is kept fixed. The first and calculable, although not necessarily dominant, 
among these corrections is of the order 0(m~ 2 ) and is given in the lowest order in a s by 
eq.(|3]). The presence of the distance scale m~ l in the 'penguin' type decays of the b hadrons 
may have further implication for the photon radiative decays as well as for the nonleptonic 
decays of these hadrons. 

I am thankful to Mikhail Shifman for useful discussions and to Mark Wise for pointing 
out a numerical error in an earlier version of this paper. This work is supported, in part, by 
the DOE grant DE-AC02-83ER40105. 
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Figure caption. 

Figure 1. A graph for the 'gluon-photon penguin'. The heavy dot represents the four- 
fermion weak interaction. At small momentum k of the gluon the scale for the expansion in 
powers of k is set by m c . 
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Figure 1 
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